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Abstract 

The polarization of Af, baryon has been measured in ALEPH, OPAL and DELPHI experiments. 
A significant loss on the transfer of the b quark polarization to the A;, baryon polarization has been 
noticed. This implies that the hadronization effects can not be neglected. Therefore, we may make 
use of the polarization measurements to look for suitable model for the A^ distribution function. 
To investigate nonperturbative effects implied in the A^ polarization, we construct four models 
based on a perturbative QCD factorization formula. The models are the quark model (QM), the 
modified quark model (MQM), the parton model (PM), and the modified parton model (MPM). 
The modified models mean the models having transverse degrees of freedom and the associated 
soft radiative corrections having been resummed. The quark and parton models can not describe 
all experiments in the same time. On the other hand, the modified models can have the power to 
explain all data in the same formalism. 

PACS numbers: 13.20Hc, 12.38Cy, 13.25Hw, 12.38-t 
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I. INTRODUCTION 



The polarization of bottom baryons A;,'s has been measured by ALEPH |]T[, OPAL 
and DELPHI The ALEPH data showed that the A;, polarization has value 

P = -0.23l|]:2o(stat.) ± 0.08(syst.). The OPAL data indicated the polarization P = 
-0.56l[!:^5j(stat.) ± 0.09(syst.). The DELPHI experiment gave P = -0.49^[}Jg(stat.) ± 
0.17(syst.). Although these three measurements are compatible with each other, the A^ 
polarization still has a wide range of value from +0.01 to —0.79. To improve the situation, 
it is better to find out a sensitive measurable quantity on the polarization. However, this is 
very difficult before we can have a more qualitative understanding on the spin properties of 
Afy baryons. This paper is intended to understand the behind mechanisms by constructing 
physical models based on perturbative QCD formalism. 

Measurement of a large longitudinal polarization of the A^ may indicate the polarization 
of primary b quark produced from a decay. The b quarks produced in the reaction 
e^e^ ^ ^ bb are highly polarized with polarization P = —0.94 [^, ^, The corrections 
from hard gluon emissions and mass effects can change the polarization of the final state b 
quarks by only 3% [0, ^. The b quark can fragment into mesons and baryons. The decays of 
b mesons into spin zero pseudoscalar states do not retain any polarization information. The 
hadronization to b baryons might preserve a large fraction of the initial b quark polarization. 
In the heavy quark mass limit, the spin degrees of freedom of b quark are decoupled from a 
spin-zero light diquark. The initial polarization of b quark can therefore be preserved until 
the Af, decays. The higher mass b baryon states can decay into the A;, baryon but transfer 
little spin degrees of freedom. These effects have been estimated from different scenarios as 
about 30%. This leads to that the final A^ polarization could be P = —0.6 ~ 0.70 [10 . 



The ALEPH Collaboration measured the A;, polarization by employing the method sug- 



gested by Bovicini and Randall |jTT|. In the ratio y(P) = {E^iV)) / {EpiV)) with {Ei) and 
{Ey) the averaged lepton and antineutrino energies in the laboratory reference frame, the 
fragmentation effects are largely cancelled out. Also, the spectra of the electrons and anti- 
neutrinos produced from the inclusive semileptonic decays of polarized A^ baryons are harder 
relative to the spectra of unpolarized decays. The ALEPH Collaboration proposed to mea- 
sure the ratio R(P) = ?/(P)/?/(0), which is a Lorentz invariant quantity. The A;, polarization 
is then extracted from a comparison between the measured ratio R^^eph _ -^ -^2 + Q.IO and 
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a Monte Carlo simulation ratio R'^^(P) with varying P. Because the A;, polarization is best 
defined in the rest frame of the Af, baryons, one can rewrite the ratio R in terms of the 
averaged variables in the rest frame to determine the polarization 

^ (e;)(e;)(i-r) 

{El){PI)R-{El){P*) ' ^ > 

where the star variables denote the averaged quantities in the rest frame and are calculated 
with P = — 1. The above equation can be investigated by different theoretical models. For 
example, if we apply the free quark model to calculate the star variables, we can obtain 
the polarization P = —0.23, which is closed to the ALEPH's result [|^]. The situation 
will become interesting as we apply the same model for the DELPHI experiment. The 
DELPHI experiment measured the same ratio Rdelphi ^ 1.21+°;J^ and obtained Pd^lphi ^ 
-0.49lo;3o. In the same way, it is easy to check that substituting DELPHI'S ratio Rdi^lphi 
into Eq. (|lD can derive a different value P = —0.38 in free quark model. On the other hand, 
in the same model, if we employ the OPAL's polarization p^^^^ = —0.56 into the ALEPH's 
and DELPHI'S Monte Carlo simulation ratios R^'^, we can extract the central value of 
corresponding R as Ri = 1.30 and R2 = L27, respectively. This seems to imply that there 
requires more investigations to find a consistent picture for Af, polarization. That is we need 
to find a model which can explain the experiments self-consistently. The model dependence 
in the equation for P, such as Eq. (|l]), arises from the z variable z^^y) = (-P/(p))/ Using 
the z variables, Eq. can be recast as 

It will become clear in the following sections that different models would give different values 
of ratio Zi but almost the same Zp due to the characteristics of the lepton and anti-neutrino 
spectra. 

In order to explore the mechanisms controlling the spin properties of polarized A;, baryons, 
we shall investigate four models. They are the (free) quark model (QM), the modified quark 
model (MQM), the parton model (PM), and the modified parton model (MPM). The parton 
model describes the probability of finding the b quark carrying a momentum fraction of the 
momentum of the A;, baryon by a parton distribution function. The quark model assumes 
that the A;, baryon contains only one b quark and two light quarks and the corresponding 
parton distribution function is just a delta function of the momentum fraction. This means 



that the b quark carries almost all the A;, baryon momentum. The modified quark model 
and the modified parton model mean that the quark model and the parton model contain 
additional Sudakov form factor and transverse momentum. The Sudakov form factor arises 
from a resummation over radiative corrections of soft gluons and have the effects to enhance 
the perturbative QCD contributions. 

We emphasize the importance of transverse degrees of freedom of partons inside a A;, 
baryon in our analysis. First, the transverse momenta regularize the divergences when the 
outgoing c quark in the process b — > cW is approaching the end point. Second, the trans- 
verse momenta also enhance the contributions from the spin vector along the polarization 
direction. For completeness, wc also introduce the intrinsic transverse momentum for the 
distribution function. We assume that the form of the intrinsic transverse momentum part 
of the parton distribution function can be parameterized as exp [—tM'^b'^] with an impact 
parameter b, which is the conjugate variable of the transverse momentum. The other factors 
are the A;, baryon mass M and a dimensionless parameter t. The impact parameter b will 
be integrated out in our PQCD formahsm. The z variables are functions of t. To determine 
the parameter t we rely on the experiments. The OPAL Collaboration determined the po- 
larization by comparing the measured distribution of the ratio E^/E^ against a simulation 
of this ratio using the JETSET Monte Carlo event generator. The polarization 'P°^^^ of 
OPAL experiment and results of the ALEPH and DELPHI experiments will determine the 
range of parameter t. 

The arrangement of our paper is as follows. In the next section, we shall demonstrate the 
factorization formula for the inclusive semileptonic decay of the A?, baryon. In this formula, 
the hard scattering amplitude, describing the short distance subprocess h civ^ convolutes 
with a jet function and an universal soft function. For simplicity, we shall assume that the 
charm quark mass can be ignored. That is we shall neglect the corrections like m^/M^ with 
TTT-c the charm quark mass. This approximation is less than 10% and is safe as compared with 
the accuracy of the experiments. However, it requires to consider the coUinear divergences 
due to our ignorance of the charm quark mass. The jet function is then necessary for 
absorbing the collinear divergences. The universal soft function involves the 6 quark matrix 
element. The matrix element contains a large scale factor, the b quark mass M^. To have 
a well established matrix element, we need to employ heavy quark effective theory (HQET) 
to scale out this large scale. We also need to separate the leading order matrix elements in 
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1/M(, expansion from the higher order ones. We shall develop a description for separating 
the leading order from the higher order mass corrections. This description is equivalent to 
the OPE approach. In Section 3, we shall construct four models based on the factorization 
formula. Section 4 gives the numerical result. The conclusion is given in Section 5. 



II. FACTORIZATION FORMULA 



We shall investigate the quadruple differential decay rate for polarized baryons, — > 

^ \Vcb?Gl 

dEidq^dqod COS Oe 256tt'^M ' ^ ^ 
where M denotes the mass of the A;, baryon, L'^" represents the leptonic tensor 

L'^'^ = Trf^^r^ ^pF'^] , (4) 

and Wf^u means the hadronic tensor 



W,^ = (27r)35(^)(PA, - q - Pxj5^(A.|crt6|X,)(X,|cr.6|A,)-^-^ (5) 



XA6|ci;&|A,)(A,|ci>|Ab)- ' 

X 

where denotes the V — A operator 7/^(1 — 75), \Xc) means the hadronic states containing 
a charm quark and q is total momentum carried by the lepton and anti-neutrino. We 
choose the normalization for A^ state jA^) as {AbiPki^ S)\MPAb^ 3)) ^ (27r)^(2P^J5^(-S - 
S')6^{P\^ — -Pa J. The kinematical variables Ei, q, go and cos^^ are expressed as follows. We 
choose the A5 baryon rest frame such that the initial A5 baryon momentum, Pa,,, and the 
final state lepton and anti- neutrino momenta, pi and pp, can be defined as 

Pa, = -^(l,l,0),p, = (p+,0,O),Pp = {pt,P^,Pu±) . (6) 

The variables E^, q, and go are related to p'l-,pt-,P9 as E^ ~ p^/-\/2, q^ — 2pfpp, and 
Qo = (pi + Pp + Pp)/a/2- We let Pf, = Pa, — / represent the b quark momentum whose 
square is set as ~ with Mb the b quark mass. The momentum / of the light degree of 
freedom of the A;, baryon can have a large plus component and small transverse components 
1_L. The final state charm quark momentum is equal to Pc = Pa, — I — q. The angle 9^ is 
defined as the angle between the third component of p^ and that of the b quark polarization 
vector, Sb- The differential decay rate can also be rewritten in terms of Ep, y, i/q and 9^ 



with Ep the anti-neutrino energy and corresponding angle Op. Because the right hand side 
of Eq. (|^) is independent of which parameterization for the leptonic variables, we shall use 
both parameterizations for the differential decay rate. 

It is convenient to use the scaled variables x^(p) = 2i?£(p)/M, y = q'^/M^, and yo = 2qo/M. 
The integration regions for y and yo are 

y 

< Xi(p) < 1, <y < h Xi(p) <yo<l + y- (7) 

xe{p) 

Note that we have chosen M as the scale factor. For simplicity, we have chosen rric = and 
left rric 7^ to other work. This approximation is safe as compared with the accuracy of 
available experiments. The leading corrections to this approximation is of order 0(m^/M^) 
being less than 10%. 

By optical theorem, the hadronic tensor W'^'^ can be related to the forward matrix element 
j^/iu tiirough the formula 

= JJ^. (8) 

TT 

The lowest order of T^*^ is defined as 

T^^(PA„g,^) = -^| d'ye"^-y{A,iPA„S)\T[f^iO),r{y)]\A,iP^„S)) (9) 

with J'^ = g7^(l — 75)6 the V — A current. The forward matrix element can be expressed 
in the momentum space 

T'^^IPa,, q, S) = -ij ^Tr[5^^(n - g)T(PA„ 5, P,)], (10) 

where the trace is taken over the fermion indices and color indices, the hard function S^'^{Ph— 
q) describes the short distance decay subprocess, h Wc, and the soft function T{P/^^, S, Ph) 
denotes the long distance matrix element 

T(PA„5,n) = J ciV^^-^(A.(PA„5)|6(0)6(y)|A5(PA„5)). (11) 

Because we are only interested in the leading contributions in this note, we are required 
to separate the leading contributions from subleading contributions. To specify the leading 
contributions, we also need to consider correction terms. The correction contributions may 
come from radiative correction terms like a" and power correction terms like 1/M™ for 
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n,m > 1. Among the radiative corrections, the contributions from soft gluons have loga- 
rithms as log(m^/(5^) will become dominate at the end points, at which the final state quark 
is approaching on-shell. As discussed in the Introduction, the hard gluon emissions can only 
contribution about 3%. Therefore, we shall retain the soft gluon contributions at the end 
points. To discuss the power correction terms, we need to be more careful. As investigated 
in the OPE and HQET approach, the power corrections can have two sources one from the 
short distance expansion for the forward matrix element and the other one from the heavy 
quark mass expansion for the expanded matrix elements. Here, we shall present a different 
approach in which the leading order matrix elements are in terms of nonlocal heavy quark 
currents composed of heavy quark effective fields in HQET. 

We now demonstrate this description. To start up, we express the forward matrix element 

tT>^%PA„ q,S) = J ^jAS'-'iPb - q)T{PA„ S, P,)] 

+ J^.J |!^Tr[5r (n - q, k')T-{PA,, S, n, k')] + ■■■ (12) 

by including a higher order term from triple parton matrix elements containing gluon fields 

T'^{PA^,S,P,,k') 

= Jd'yJ dSe^^^V(^^-'=')-(A,(PA,,5)|6(0)(-i^A-(z))%)|A,(PA,,5)) (13) 

with A°'{z) the gluon fields. We shall employ the hght-cone gauge — n ■ A — 0. To 
continue, it is useful to introduce the light-cone vectors p and n in the -|- and — directions, 
respectively. These two vectors satisfy properties p"^ ~ — and p • n — Pa^ • n. The A^ 
baryon momentum P\^ is then recast as 

= + n" . (14) 

For the b quark inside the A^ baryon, we parameterize its momentum P^ as 

= ^P" + ^i^ri^ + PL (15) 

p2 _ /Vf 2 

where P^ = is the on-shell part of Pj, and the momentum fraction z defined hj z = 
P^ / Pj[ = 1 — l~^/P^ . By the parameterization of P^, the b quark propagator is then 
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expressed as 



^b-Mb + ie P'^ - M2 + ie 2Pb ■ n + te 



2 ]\/r2 I o o „!„■,• 



Since we allow ^ P^^, therefore we have 2; ~ A/M with A about Aqcd and ~ M. The 
second term of Eq. (|l^ having power like ^ ~ is as large as the first term. We thus do 
not need to take into account the power correction from the coUinear part of the b quark 
propagator. This demonstration can be better understood from the fact that the b quark 
is almost on shell and has a large quark mass. Therefore, there is no collinear divergences 
associated with the b quark. The remaining work is to separate the leading terms of the 
hard function 5*^*^ from the higher order terms. The hard function S'^'^ is a function of I. 
We can make Taylor expansion for S'^'^ with respect to l^. This is because the momentum 
/ has a large plus component = with ^ = 1 — z. By performing Taylor expansions for 
and S^^l.k') around S'"'{ip) and S^" {^p.^'p), we then obtain 

S^^{1) = S^%ip) + ^U,{l-ipr + -- - , 
The Ward identity ensures the following equation to hold 



(ep) = 5r(ep,ep). (19) 

The contributions from S^^ terms are power suppressed than S^'^{C,p) by at least 0{1/M^). 
The effects of the second terms in Eq. ( |18D are to replace the gluonic field operators in the 
second term in Eq. ([I^ ) by covariant derivative operators. Let's explain this. By adding the 
second terms from Eqs. (0) and (p!2D , respectively, we have 

/ +/ (^/ ^^r(eP,eV)T"(/,A;') . (20) 

In light cone gauge n ■ A = 0, we can rewrite the above equation as 

/^/ ^^r(eP,e'pXrT(/)(27r)^5^(A;') +r°'(/,A;')] (21) 

where the projection tensor w"' = g^' ~ has been employed. Using the identity 

{27i)''6\k') = j dSe^^'-^y-'^ , (22) 



the bracket in Eq. ( pi]) can be recast as 

dy' [ d^ze^'V'='-(^-^)(A,(PA„5)|6(0)D"(z)%)|A,(PA„5)) (23) 



with = id"" - igA". 

The b quark field in the leading matrix element T still contains a large phase 
exp {—iMbV ■ x) with v the Af, velocity. This is unsuitable to define a matrix element at 
low energies. To solve this, we can employ the heavy quark effective theory (HQET). In 
HQET, We can rescale the b quark field, b{x), into by{x) = exp {iMbV-x)^'h{x). The 
rescaled b^ field is a small fluctuation quantity of coordinate, since the remaining scale in 
its phase factor is only about Aqcd scale. In HQET, Pb is parameterized as Pb = MbV + /c, 
with k the residual momentum. The rescaled b quark field, by{x), carries the residual mo- 
mentum k and has a small effective mass A, with A = M — Mb. Under the heavy quark 
mass expansion 

b,{x) = K{x) + 0{^) + --- , (24) 

the matrix element T in terms of 6^ can be expanded as 

^ = ^o + 0(^) + --- . (25) 

The To is in terms of an universal effective heavy field h^, which is defined as the 6„ field 
in the infinite mass limit Mb — > oo. The missing of 0{j^) term is due to the equation of 
motion. The expression for Tq is easily written down as 

To= f d'ye''-y{Ab{v,S)\Kmv{y)\Hv,S)) . (26) 



Note that we have replaced the hadronic state vector |A;,(PAi,) S)) by its equivalent represen- 
tation \Ab{v,S)). The normalization of \Ab{v,S)) is large than the usual normalization by 
a factor M^. We shall skip the derivation on how to derive the above equation. To derive 
leading order contributions, we still need to extract the leading spin structure of Tq. This 
can be achieved by means of Fierz identity. As a result, the leading order forward matrix 
element T^*^ takes the form 

T'^^(PA„g,^) ^ ~^| 0^^!^[S>^-'{k = Cp,q)mTo{PA„S,k = ^p)^^] 

-[S^^k = Cp, q) ^,75][To(Pa„ 5, k = iP)^^]\ + O(^) (,27) 



where we have inserted the Fierz identity 

'^l5)im{l5lfj)jn^ (28) 

where i,j,m,n denote the fermion indices and the dots represent the other gamma matrix 
would result in higher order terms. 

We now briefly describe how to derive the factorization formula for the inclusive semilep- 
tonic decay A5 Xgiu. The details about the derivation of the following factorization 



formula can be found in |12|. We shall only demonstrate the main ideas and not try to give 
a repeated proof. The formula for the quadruple differential decay rate can be expressed as 
a convolution integral over the soft function 5*, the jet function J and a hard function H 

1 d^T M2 /-^-^ , f ,2, 

dz / d I 



r(°) dxdydyod cos 9 2 

xS{z,\^,fi)J{z,p-,l^,fi)H{z,P~,fi) , (29) 



where x = Xe{xp ) andO = 6^(6^,) and T^o) = ^\Vgh\^M^ The scale /i is introduced as a 
renormalization and factorization scale. The transverse momentum 1 1 has been reintroduced 



for regularization of the end point singularities |T2[. The end point singularities arise from 
the end point region x —>■ 1 and y, y^ —>■ 0. The charm quark (assumed as massless) 
has a large minus component = (1 — y/x)M/\/2 and a small plus component P^ = 
(1 — yo — ylx)Mj\/2. This implies there is a very small invariant P'^ = M^(l — yo + y), 
which leads to an on-shell jet subprocess. The 1_l integrals can be finished only when we 
know the exact dependence of the jet function on 1^. But the jet function is nonperturbative 
and cannot be determined theoretically, so far. Fortunately, this difficulty for integration 
over l± can be removed by means of a Fourier transformation for the jet function into its 
impact space representation as 

J(^,p-,U,/i) = J i!^J>,p-,b,/.)e^'-^ . (30) 

The \± integrals then decouple from the jet function and the remaining factor e*'^ "^ is then 
associated with the soft function. The factorization formula Eq. ( p9l) can also be applied to 
the case with loop corrections. With the Fourier transformation for I^, the Feynman rule 
for the radiative gluon cross over the final state cut should be modified with an extra phase 
factor e*'^ "^. The upper and lower limits of z are chosen as Zmax = 1 and z^^^^ = x. The 
lower limit z^am = x is from the jet function. The upper limit z^i,^ = 1 is chosen to fill 
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the kinematical gap between Mi, and M. The Fourier transformation of Eq. 
impact b space then takes the form 



91) into the 



r(°) dxdydyod cos 9i 



M2 



dz 



d^h 



(27r)2 

X S {z, b, iJ,)J{z, P~, b, (z, P^, /i) 



(31) 



To deal with the colhnear and soft divergences resulting from the radiative corrections for 
massless parton inside the jet, the resummation technique is necessary and these divergences 
could be resummed into a Sudakov form factor |T^. The jet function is then re-expressed 
into the form 



J{z,Pc ,b,fi)=exp [-2s{P^ ,b)]J{z,b,fi), 



(32) 



where exp [—2s{P^ , b)] is the Sudakov form factor. The RG invariant Sudakov exponent has 
the expression up to one loop accuracy 

AW 



In (26) 



2/?i 

1 In (2g) + 1 



(q-b) 



ln(- 



ln( 



with the variables 



-^[\n\2q)-\n\2b)] 



(33) 



(34) 



We choose the QCD scale A = Aqcd to have the value 0.2 GeV in the numerical analysis in 
section 4. The other factors are defined as 



33 - 2n 



f 



12 ' 
153 - 19n/ 



24 



3 
4 



a(2) 67 7r2 10 8^ , ,e\ 

= Uf + -/?iln(— ) . 

9 3 27 ^ 3^^ 2 ^ 



(35) 
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The scale invariance of the differential decay rate in Eq. (^T]) and the Sudakov form factor 
in Eq. (p^) requires the functions J, S, and H to obey the following RG equations: 



-D Jib, fx) 
M(P,-,/i) 



-2jgJ{b,fx) , 
-lsS{b,fi) , 



(36) 



with 



d d 



(37) 



7g = —as/n is the quark anomalous dimension in axial gauge, and 75 = —{as/7i)CF is 
the anomalous dimension of S. After solving Eq. (p6|), we obtain the evolution of all the 



convolution factors in Eq. (|3T| ) 



Pq ,b,fi) = exp 
S'(2, b, /i) = exp 

H{z,P~,fM) = exp 



-2.(p-,6)-2/ ^7>.(/^)) 
i/fe 

/(^,&, 1/&) 



J(z,6,l/6) 



—IsiasifJ')) 

i/b 



[27g(a.(/^)) +75(tts(/^))] 



H{Z,P-,P,'). (38) 



In the above solutions, we set the 1/6 as an IR cut-off for single logarithm evolution. 

For the initial soft function f{z,b,l/b), we shall keep the intrinsic b dependence by 
f{z,b,l/b) ~ f{z,b). The b dependence in f{z,b) can support us a way to explore its 
effect in determining the polarization. We assume f{z, b) to have the form 



/(z,6)=/(z)e-^('') 



and take an ansaze for parameterizing exp [— S(6)] as 



(39) 



(40) 



with an unknown parameter t. To avoid double counting for the contributions from trans- 
verse degrees of freedom, we need some modifications for the factorization formula. For the 
end point regime where the Sudakov suppression dominates, we employ the approximation 



f{z,b)=f{z) 
12 



(41) 



, while for other regimes which are not under the control of the Sudakov suppression, we 
take into account the intrinsic b dependence of f{z,b) 

f{z,b)=e-"^'''f{z) . (42) 

We make further approximations such that J{z,b,l/b) = J^^\z,b), and H(z, , P^) = 
H^'\Z,P-). 

Combining the above results, we arrive at factorization formula as 
1 d^T 



r(o) cixdydyod cos 9 



f dz r^jW(z,6)ij(°)(^,p-)e-^(^^"''') 
Jx Jo 47r 

f(z) for X in end point regimes , 

exp[—tM'^b^]f{z) for x in other regimes , 



S{P-, b) = 2S{P-, b) - r ^[27,(/i) + ism ■ (44) 

Ji/b 



where 

h/h 

The parameter t will be determined by experiment. From practical calculations, we find that 
the above difference between the distribution function with and without intrinsic transverse 
momentum contributions is very small. Therefore, we shall include the factor exp[— tM^6^] 
for entire range of x in the numerical analysis. 

Let's now discuss how to parameterize To(fc) defined in Eq. ([TT|). As discussed in previous 
paragraphs that, at leading order, To(/c) is expanded in the form 

T,{k) = -nj/fiz) - -^S, ■ njfy.giz) . (45) 

The unpolarized and polarized distribution functions, f{z) and g{z), are defined as 

= / ^e^^'-'^'-{v,S\K{0)r/Ki^)\v,S) (46) 
J 2n F ■ n 

and 

g{z) = [ ^e^(^-^^'-{v,S\K{0)r/lM^)\v,S) . (47) 
J Ztt F ■ n 

It is easy to show that f{z) and g{z), in the heavy quark limit, share a common matrix 
element which could be described by an universal distribution function, fA^{z). This just 
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reflects the heavy quark spin symmetry. We adopt the distribution function proposed in |T2 
in the form 

The parameters A^, a and e are fixed by first three moments of fivSz) 



fA,iz)dz = 1 , 



b 
1 



dz{l - z)fAS^) = A/M + 0(A^cd/M' 

1 a2 



2 dz{l - zfh,{z) = ^2+l^b + 0{AlcJM') , (49) 
where A = M — and is to parameterize the matrix elemnet 

By substituting the inputs 

M = 5.641GeV, M,, = 4.776GeV , i^b = 0.012 ± 0.0026 , (51) 
into Eq. (H), we determine the parameters A^, a and e to be 

N = 0.10615 , a = 1 , e = 0.00413 . (52) 
For simplicity we shall omit the subscript of fK^{z) in the following text. 

III. DIFFERENTIAL DECAY RATES 

In this section, we construct four models based on the factorization formula Eq. (^). The 
models are the quark model (QM), the modified quark model (MQM), the parton model 
(PM), and the modified parton model (MPM). The charged lepton and anti- neutrino spectra 
for the decay Af, Xgiu in the quark model are expressed as 

1 ^^r^^ ^ I f [(3 - 2xe) - Pcos9e{l - 2x,)] for i 
T(<^Uxdcos9 I ^2;p(l-Xp)(l-Pcos^p) for z/ 

where P and cos 6'^(p) denote the polarization and cosine of the angle 6'£(p) between the third 
components of the lepton (anti-neutrino) momentum and the A;, spin vector. 
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By taking into account Sudakov suppression from the resummation of large radiative 
corrections, and substituting f{z,b) = 5{1 — 2;) exp [— tM^6^], H^^^ = (xi — y)[{yo — xe) — 
P cos Oi{yQ — X£ — 2y/xi)] and the Fourier transform of J*^°) = 6{P'^) with P'^ = M^(l — ?/o + 
y ~ Pl/^b) i^^*^ -^1- (EH)' derive the lepton spectrum in modified quark model . The 
spectrum is, after integrating Eq. (E3|) over z and yo, described by 



r(°) dx id cos 9, 

X 



/o Jo 

{{l + y- Xi) - P cos0f(l + y-Xi- 2^))Ji{r,Mh) 

Xi 

' ^ riJ2[riMh)-r]^J'i[riMh)){\-VcosQi 



Mb 



(54) 



where P^^ = {l—y/xe)M/y/2, rj = ^y {x£ — y){l /x£ — 1) and Ji,J2,<^3 are the Bessel functions 
of order 1, 2 and 3, respectively. Note that we have made an approximation by substituting 
exp [—tA'Pb'^] for the end point regimes. We also need the neutrino spectrum in modified 
quark model as 

^ ^'^m^ = M r dy f"^ dhe^-'^^^^^^'X-'^P^^^x^il - x.WMv'Mb) 



r(°) dxpdcosOp 



^0 



x(l-Pcosep) (55) 



with T]' = -y){l/xp - 1). 

The charged lepton spectrum in parton model is obtained by adopting H^^^ = [xg — 

y)[{yo-xi, - (1 - z)y/xi>) -Fcos9e{yo- Xi - i^ + z)y/xi)] and = M'^[l - y^ + y - (1 - 
z){l — y/xe)]. After integration over yo, we then derive 



1 d^T 



PM 



r(°) dxid cos 6 f 



/ dy dzf {z){xi - y)[{y + z - Xi) -F cos 9i{y + z - Xi - 2z—)] . (56) 
Jo Jxe 

In the same way, the neutrino spectrum can be written down 

I dy I dzf{z)xp{z — Xp){l — P cos 6p) . (57) 

Jo Jxa 



1 r^, ri 



PM 



PW dxpdcosOp 

The charged lepton spectra in the modified parton model takes into account large per- 
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turbative corrections and nonperturbative intrinsic contributions with the expression as 

^ : M r dy t dz t"^ d6e[-*"™^^''\-^(^-'') f{z){x, - y)v 



r^°'^ dxed cos 9 i jq j^^ jq 



X 



{{z + y- xA - P cosOeiz + y - xe - 2z—))Ji(r]Mb) 

Xi 

' ^ -r]J2{vMb) - Tf J^{r]Mh)){l - P cosOi) 



(58) 



W6 

with Tj = {x — y){z/x£ — 1). The neutrino spectrum in modified parton model is also 
easily derived as 

i a i MPM , , / ,1,1 



dxpd COS 9p 



a'li 



M/ dy dz rf6e-^(^^"'^)e-*"™^'^V(^)a;p(^-Xp)VJi(VM6) 



x(l-Pcos^p). (59) 
with r]' = a/(xp -y){z/xp - 1). 

IV. NUMERICAL RESULT 

The Afc's produced in ALEPH, DELPHI and OPAL experiments are highly boosted in 
the laboratory frame. For the relativistic A^'s, the forward-backward asymmetry of a decay 
product can be directly expressed in terms of a shift in the average value of its energy. The 
charged lepton also carried a residual sensitivity to the Ab polarization. Because neither the 
Ab four-momentum nor the lepton four-momentum can be fully reconstructed in the exper- 
iments, the ALEPH and DELPHI experiments proposed to measure the A;, polarization, P, 



through the variable y suggested in |11 



<Ee> 

However, there still exist many uncertainties suffered from experimental procedures on ex- 
tracting the energy spectra. It requires normalizing the measured y with an unpolarized 
simulated ?/*^'"(0). Therefore, the experimentally measured quantity is the ratio 

y{P) 

ALEPH and DELPHI determined the polarization by comparing the measured value of 
ratio R from the Monte Carlo simulation R'^'^(P) with varying P. The experimental results 
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are R = 1.12 ± 0.10 and P = -0.23+E}:^^(stat.) for ALEPH and R = 1.2ll°:J^ and P = 
—OA9'to^l{stat.) for DELPHI, respectively. Theoretically, the A;, polarization can be best 
defined in the rest frame. It is instructive to rewrite y in terms of average variables in rest 
frame 

^ {Ei)+p{pnp = -i)) . . 

^ (^;:) + p(^.*(p = -i)) ' ^ ^ 

where the star average variables are evaluated with P = — 1. The average variables can be 
calculated from the formula 

f f „ dxd cos 9 , , 

/ \ J J dxd cos (P.'i\ 

by employing different models for the differential decay rate. It is much simplified in cal- 
culations of these averages, if the charged lepton and anti-neutrino average quantities are 
evaluated by their corresponding models for the differential decay rate. From these relations 
we can determine P in terms of R as 

^ {E!){E:)il-R) 

(E;)(p^)R-(Ei)(Pr) ■ ^ ^ 

We first compare the difference between the experimentally determined polarization P^^^ 
and the theoretically evaluated polarization P^^ in the four models QM, PM and MQM, 
MPM with parameter t^^^ = t^™ = 0. The result is shown in Table. I. We can see that 
the theoretical polarizations are close to the ALEPH polarization P^leph _ _q 23 but have 
a large deviation from the DELPHI polarization p^elphi _ _q 49 Among different model 
evaluations with one R, their differences are very small. This implies that nonperturbative 
effects from distribution function over longitudinal momentum fraction and perturbative 
effects from Sudakov suppression are not important in determining the polarization. 

We now turn on the parameters t^Q^ and t^^^ to find out their values from experiments. 
It is interesting to note that the ratio zi = {PD j {ED is model dependent but the ratio 
Zp = (Pj?) / (Ep) is almost the same for all models. Using these two z variables, we can 
rewrite Eq. (|^) as 

1-R 



ZpH — Z£ 

Since 2;p ~ 1/3 for all models, we can further simplify the above equation into 

p _ 3(1 - R) 
R — 3zi 
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(65) 



(66) 



We show the behaviors of P with respect to for Ral^ph ^ ^ y^. ± 0.10 (ALEPH) and 
^DELPHI ^ 1.21+°;^^ (DELPHI) in Figs. 1 and 2, respectively. By applying the experimental 
bounds for P, we can extract from Fig. 1 the zi range of — oo < zi < —0.105 for ALEPH 
and from Fig. 2 the range of —1.75 < < —0.02 for DELPHI. Theoretically, the Zn range 
would be model dependent. Considering MQM and MPM and plotting the z^ — t relation 
in Fig. 3, we can find that the maximum of Z(, can not be larger than —0.05 with t ~ 0.3. 
This is because the suppressions from the contributions of intrinsic transverse momentum 
are modeled by parameter t. By varying the value of we can easily change z^. However, 
the fluctuations from the Bessel function in the differential decay rates would prevent the 
suppression of t from becoming large. In the end, there exist a maximum bound for z^. We 
notice that there is also a lower bound for zi as zn > —0.18 with t ~ 2. We also hope that 
t should be less than unity and close to zero to make the perturbative calculation reliable. 
Combining the above considerations, we can obtain the range for Z£ as —0.12 < ze < —0.105. 
The lower bound of ze comes from t — 0. From Fig. 3, we can employ the z^ range to derive 
the t range. The ze — t relation is a bounce with maximum at t ~ 0.3 for both MQM and 
MPM. The Zi range implies that < t < 0.1. We have taken the convention that t is located 
left to the maximum. This convention for choosing the value of t can be further checked 
after we discuss the OPAL experiment. The Zi — t plots for MQM and MPM are not much 
different. 

We now discuss the possible constraint over ze can be obtained from the OPAL experi- 
ment. The OPAL Collaboration employed a comparison between the measured yS — Ep/Ee 
and the Monte Carlo simulated yf^ to determine the polarization P = —0.56to^l{stat.). 
Applying the OPAL P to DELPHI and ALEPH experiments, we obtain —0.6 < ze < —0.1 
for DELPHI and -0.55 < z^ < -0.105 for ALEPH. The bound of ze is -0.55 < ze < -0.105. 
Looking at Fig. 2, one can see that the OPAL experiment give bound on Zj^ in consistency 
with previous investigation. In summary, the ALEPH, OPAL and DELPHI experiments 
imply the range value of < i < 0.1 and the corresponding range of —0.12 < z^ < —0.105. 

As a consistent check, we can write R in terms of P and z^ as 

R="^ylf. (67) 
(3 + P) ^ ^ 

By this equation, we can parameterize the Monte Carlo simulation ratios it!'^*^(P)'s of ALEPH 
and DELPHI. We flnd that the value of z^ ~ —0.075 can be used for both experiments 



18 



in a good approximation within 5 — 10%. In Fig. 4, we compare the R — P plots for 
Zi ~ —0.075 and —0.12 < ze < —0.105. The experimental bounds for ratio R can give 
constraints over P. The combinnation of ALEPH and DELPHI experiments gives bound 
of P as -0.79 < P < -0.05, while the theoretical bounds are -0.73 < P < -0.05. The 
theoretical bounds for P being smaller than the experimental ones can be easily understood 
from the maximum bound of Z£ from theory is smaller than the Z£ employed in the Monte 
Carlo simulation performed by experiments. The difference between theory and experiment 
can be compensated by including higher order corrections, such as the mass corrections, etc. 

V. CONCLUSION 

In this paper we have constructed four models based on the PQCD factorization formula 
for Ab — > XJu. We used these models to investigate the physics implied by the ALEPH, 
OPAL and DELPHI experiments. We found that these experiments can be understood 
from theoretical models, the modified quark model and modified parton model. These two 
models contains intrinsic transverse momenta for partons, which are nonperturbative and 
parameterized by an exponential form with a parameter t. The parameter t relates to the 
variable Zi = (P/(P = — with (P/(P = — 1)) and (E^) the average momentum 
and energy of charge lepton in the rest frame of A;, baryon. Wc found that the ratio 
R = y{P)/y{0) can be approximately expressed in terms of P and ze. Using experimental 
results, we then determined the ranges of z^ and t. 
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Table. 1 The values of the A;, polarization are predicted from the quark model, the parton 
model, the modified quark model and the modified parton model by employing the ALEPH 
and DELPHI experiments. The ALEPH and DELPHI experimental results are also shown 
for comparison. 



Pqm PpM Pmqm Pmpm 


Pexp 


R 


Experiment 


n OO+0.19 Qo+0,19 n Q 1+0.20 n Q 1+0.20 
'-'•-'^-0.17 ^-'■-'-'-().i7 ^-'■--^-0.17 -'-0.17 


*-'■-'-'-(). 20 


1.12 ±0.1 


ALEPH 


n qo+0.24 r, qo+0.24 r, qQ+0.25 n qQ+0-25 
U.OO_Q 24 U.OO_Q 24 U.Oy_Q 24 u.oy_Q 24 


_n 4Q+0-32 
"-"•^^-0.30 


1 91+0.16 
^-•^-■--0.14 


DELPHI 
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FIG. 1: Plot of P vs. zi is shown by employing the ALEPH ratio R = 1.12ib0.10. The experimental 
polarization P = — 0.23^q2o is also shown to indicate the allowed range of z^. 
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FIG. 3: Plot of Z£ vs t. The modified quark model (solid line) and modified parton model (dashed 
line) are shown. 
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FIG. 4: Plot of R vs. P. The ALEPH and DELPHI Monte Carlo simulations (solid line) and the 
theoretical prediction (band line) are shown. 
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